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ABSTRACT 


We discuss the detectability of gravitationally bounded pairs of gas-giant 
planets (which we call “binary planets”) in extrasolar planetary systems that are 
formed through orbital instability followed by planet-planet dynamical tides dur¬ 
ing their close encounters, based on the results of N-body simulations by Ochiai, 
Nagasawa and Ida (Paper I). Paper I showed that the formation probability of a 
binary is as much as ~ 10% for three giant planet systems that undergo orbital 
instability, and after post-capture long-term tidal evolution, the typical binary 
separation is 3-5 times the sum of physical radii of the planets. The binary 
planets are stable during main sequence lifetime of solar-type stars, if the stellar¬ 
centric semimajor axis of the binary is larger than 0.3 AU. We show that detecting 
modulations of transit light curves is the most promising observational method to 
detect binary planets. Since the likely binary separations are comparable to the 
stellar diameter, the shape of the transit light curve is different from transit to 
transit, depending on the phase of the binary’s orbit. The transit durations and 
depth for binary planet transits are generally longer and deeper than those for 
the single planet case. We point out that binary planets could exist among the 
known inflated gas giant planets or objects classihed as false positive detections 
at orbital radii > 0.3 AU, propose a binary planet explanation for the CoRoT 
candidate SRcOl E2 1066, and show that binary planets are likely to be present 
in, and could be detected using Kepler-quality data. 

Subject headings: occultations - techniques: photometric - planets and satellites: 
detection 
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INTRODUCTION 


i£l. 


More than 7000 transit-like light curves have been obtained by Kepler observation^ 
Among them, more than 3000 objects are identihed as planetary candidates and more than 
2000 objects are false positives. Over 800 of them show transit depth comparable to that 
caused by gas giant planets. In this paper, we point out that some of the known inflated 
gas giant planets and objects classified as false positive detections could be gravitationally 
bounded pairs of gas-giant planets, which we call “binary planets.” 

Orbital stability of exomoons around gas-giant planets hav e been studied, and these 


moon s are stable except for those around close-in gas giants (e.g.. 


Namouni 


20 mi 


Gong et ah 


2013 1. “Habita 
discussed (e.g.. 


jility” of the moons around 


Williams et ah 


19971: 


Heller 

2012; 

Heller & Barnes 

20131 


201311 . In addition, the 


observational detectabilit y of exomoons has been discussed and many detection methods 
have been proposed 


Sato&Asada 


20091 


(e.g., 

Sartore 

tti & Schneidei 

1999; 

Simon et al. 

2007; 

KinninG 

2009a b 

Kaltenegger 

2010; 

Kinning 

2012; 

Zhuang et al. 

2012 

Lewis 

2013 

)• 


Gas giants that are distant from their host star may co mmonly have moons, 


because 


regular moons are fo r med in circumplan etary disks (e.g.. 


Ganup fc Ward 


2006 


Sasaki et ah 


Mosqueira fc Estrada 


2003alJb : 


20101 ) an d the formation of the se disks is a part of 


Ganup fc WardI (1200611 proposed that the 


gas-giant planet formation. However, since 
maximum mass of a moon may be < 10“^Mp where Mp is the planet mass, detection of 
these moons is not easy. 

On the other hand, it is relatively easy to distinguish bi nary gas gian t s from single gas 


giants, because the companions are as large as the primary. 


Ochiai et al. 


( 2014 hereafter 


^ http: / / exoplanet archive. ipac .caltech.edu / cgi-bin / ExoTables / nph- 


exotbls?dataset=cumulative 
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refereed to as Paper I) found through N-body simulations incorporating both planet-planet 
and planet-star tidal interactions, that the formation rate of binary planets is as much 
as ~ 10% of the systems in which orbital crossing among multiple gas giants occurs. 
Furthermore, Paper I also predicted the binary orbital separation distribution and the limit 
of stellarcentric semimajor axis of the binaries beyond which the binary orbits are stable 
during the main sequence lifetime of solar-type stars. Thereby, binary planets may have 
more chance to be detected than exomoons if these theoretical predictions are correct!^ 


The importan 


Ochiai et ah 


point is that the mechanism to form the binary planets proposed by 


(120141) is one of the natural outcomes of orbital crossing among gas giants. 


The gas giants so far discovered in extrasolar planetary systems often have eccentric (say, 
e > 0.2) orbits, except for close-in planets where the eccentricities are tidally damped. 

It is considered that orbital crossing for the case of a system containin g three gas giant s 


or more, is the most likely orig i n of gas giants in eccentric orbits fe.g.. 


Marzari fc Weidenschilin 


12002 


Chatteriee et af 


2008; 


Lin fc Ida 


Juric &: Tremaine 


I997I: 


20081) . Typical 


fates of such three planet systems are ejection of a planet, planet-planet collisions, and 
planet-star collisions. 


However, if planet-star tidal interaction is taken into account, most o 


star collisions are replaced by formation of “hot j up iters. 


Nagasawa et al. 


:he p lanet 


(120081 ) and 


Beauge &: Nesvornvl (120121) found through N-body simulations that hot j up iters are formed 


in as much as 10-30% of the system s. The discovery of retrograde hot jupiters (e.g.. 


Narita et ah! 


2009 : 


Winn et al. 


20091) strongly suggests that some fraction of hot jupiters 


were formed through the tidal capture. 


^If exomoons can be formed in a different way than 


Canup fc Ward! (120061) considered. 


exomoons, can be larger and their detection is not so difficult. 
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Podsiadlowski et al.l (120101) pointed out that incorporation of planet-planet tidal 
intera ctions replaces some of plan et-planet collisions with formation of binary planets. 


While 


Podsiadlowski et al. 


fl2010l) assumed arbitrary initial conditions of two giants in 


closely packed, nearly circular orbits, Paper I considered more appropriate initial conditions 
consisting of thr ee giants i n sepa rated orbits, and found that the formation rate is still as 


much as 10%. 


Ida et al. 


(120131) showed through a planet population synthesis simulation 


that such orbital crossing among gas giants commonly occurs in systems formed from 
relatively massive protoplanetary disks. 

Possible methods to detect binary gas-giant planets include radial velocity, transit 
light curves, transit timing variations (TTV), Rossiter-McLaughlin (RM) effect, and 
gravitational microlensing. Because binary planets may be tight binar ies, the detection of 


such b inary planets will be through deviation from a single planet £t. 


Podsiadlowski et al. 


(120101) showed that the radial velocity amplitude of the deviation is too small. Detection by 
gravitational microlensing is also possible if a background star passes between the binary 
planets, but they suggested that the possibility of such events is too low. We found that 
TTV signals are also too small. While detection by RM effect is not ruled out, a bright 
host star is required, and the number of sufficiently bright stars with transiting gas giants 
beyond 0.3 AU is limited. Therefore, transit (light curve) observation is the most promising 
method to detect these binary planets. 


In addition, binary planets may be sufficiently numerous to be detected in large 
transiting planet surveys. Assuming that the ~ 10% of planetary systems hosting eccentric 
gas giants are all a result of orbital instability, and combining this with the predicted 
~ 10% binary planet formation probability, we have that ~ 1% of systems could host a 
binary planet. Taking a typical stellarcentric semi-major axis of 0.5AU, this gives a transit 
probability of ~ 1%. Consequently, we expect that approximately 1 in 10000 stars would 
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host a transiting gas giant binary planet pair. This compares favonrably with the ~ 145000 
stars monitored by Kepler and ~ 100000 light curves produced by CoRoT. 


Although 


Podsiadlowski et ah 


fl2010l) also concluded that transit observations are the 


most promising way to detect binary planets, they did not discuss how planet binarity 


modifies transit light curves, or calculate detection proba 
calculate long-term tidal evolution of the binaries either. 


jilities. In add i tion, t hey did not 


Sato fc Asadal ( 20091) calculated 


the modulation of transit light curves, but as they considered Earth-Moon like systems, 
focusing on mutual eclipses, their arguments cannot be applied to the binary gas giants 
that we consider here. 


In this paper, based on the detailed theoretical predictions by Paper I, we calculate 
the effect of planet binarity on transit light curves and discuss the possibility of detection 
of binary planets. In section [21 we summarize the results of Paper I. In l|3l we predict light 
curves of possible binary planet systems and discuss the detectability of extrasolar binary 
planets. In we apply this work to the problem of detecting binary planets in the CoRoT 
and Kepler candidates, while is a summary. 


2. THEORETICAL PREDICTIONS 

Here we briefly summarize the results of theoretical calculations in Paper I. Paper I 
carried out two sets of calculations: 1) N-body simulations of three gas giants incorporating 
planet-planet dynamical tides (as well as planet-star tides) on timescales ~ 10^ years to 
investigate tidal capture to form binary planets and 2) numerical integration for long-term 
evolution of the binary orbits due to planet-planet and planet-star quasi-static tides during 
main sequence lifetime of solar type stars (~ 10^° years). 

Dynamical behaviors in two planet systems are qualitatively different from those with 
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three planets or more. For the case of two planets, they immediately start orbital crossing 
when their initial orbital separation is smaller than a particnlar critical valne, while close 
enconnters never happen otherwise. In systems with three planets or more, there is no solid 
stability bonndary. With modest initial orbital separations, three planet systems can start 
orbital crossing after their eccentricities are bnilt np over relatively long timescales (e.g.. 


Chambers et ah 


19961) . Because it is not easy to establish unstable orbital configurations of 


two planet systems. Paper I considered the three gas giant planet systems having equal-mass 
of Mj and radius 2i?j, where Mj and i?j are the present values of Jupiter. The results of 
the N-body simulations are summarized as follows: 


1. During close encounters, energy dissipation due to planet-planet tides often results 
in formation of a gravitationally bound pair of planets (binary planets). Tidal 
capture usually occurs in the early phase of the orbital crossing before the planets’ 
stellarcentric eccentricities have not been maximally excited. The formation rate is 
~ 10% of the runs, almost independent of initial stellarcentric semimajor axes of the 
planets. 

2. The stellarcentric semimajor axes of the binary barycenters are comparable to the 
initial locations. Their stellarcentric eccentricities are distributed in a broad range 
with median values of ~ 0.15. 

3. The binary planets are tight binaries. After the orbital circularization and long-term 
tidal evolution, the binary separations are ~ 3 — 5 times of the sum of the physical 
radii of the planets (7?tot)- 

Hereafter we use the subscripts ”0”, ”1”, and ”2” to represent the period at the tidal 
capture, tidal circularization, and the spin-orbit synchronous state just after quasi-static 
planet-planet tidal evolution, respectively. Because tidal interaction is a sensitive function 





of distance, the close encounters that lead to tidal capture are usually grazing ones. So, 
the binary orbits just after the tidal capture have pericenter distance of qu,o ~ Rtot — 2i?tot 
and binary orbital eccentricity ebi,o ~ 1- The binary separations after tidal circularization 
are given by abi,i ~ 2gbi,o ~ 2i?tot — 4i?tot due to conservation of angular momentum 
(Tobip ~ ^abi,o(l - e^io) ~ 


After the tidal trapping to form binary planets and the binary orbital circularization 
due to planet-planet dynamical tides, the binary separation expands, and the binary planet 
pair enters a spin-orbit synchronous state through quasi-static tidal evolution. The initial 
total angular momentum of the circularized binary is Li = 2 x {2/5)MpR‘^Up + /ia^j if^biq, 
where fi = Mp/2, Mp is the planetary mass, Rp is the planetary radius, cOp is the spin 
angular velocity of the individual planets, f2bi,i = y^2GMp/a^, and G is the gravitational 
constant. In the spin-orbit synchronous state after quasi-static planet-planet tidal evolution, 
the total angular momentum is L 2 = 2 x (2/5)Mpi?pffbi,2+ habi, 2 ^bi, 2 - For initial spin period 
of 10 hours, Mp = Mj and Rp = 2Rj, the orbital separation in the spin-orbit synchronous 
state is given by 


Obi ,2 


®bi,l 1 + 


5 V GMpabiq^P"^” 


®bi,l 1 + 0.4 


Obi,! 

Tor; 


- 0 . 5 ' 


( 1 ) 


where we assumed abi ,2 S> Rp. The distributions of abi,i and abi ,2 obtained by Paper I are 
showin in Fig. [1] At stellarcentric distance oq ^ aG.mii — 0.2 AU, Obig exceeds rn/S where 


rn is Hill ra dius (rn = (I 


the binary fISasaki et al 


/Ip/SM ^Y^^ag) and the perturbations of the central star destabilize 


2 OI 2 I) . At ttQ < OG.tide — 0.4 AU, planet-star quasi-static tide 


removes the binary orbital angular momentum and the binary planets collide with each 
other within main sequence phase of solar type stars (~ 10^° years). Since the planets’ gas 
envelopes fully contract in 10® years, in the quasi-static tidal evolution Rp = Rj is more 
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appropriate than i?p = 2i?j. With i?p = i?j, the critical stellarcentric semimajor axis is 
flG.tide ~ 0.3 AU, rather than 0.4 AU. Therefore, binary planets should be orbitally stable 
at oq ^ 0.3 AU. 

Since the timescale to establish the spin-orbit synchronous state is < 10® years and 
binary separation does not evolve signihcantly once this state is achieved, except in the 
case of the orbital destabilizations, observed binary planets should be in this spin-orbit 
synchronous state, that is, binary separations of 3i?tot^5i?tot, which is comparable to the 
stellar diameter. 

In addition, these binary systems should be stable when oq ^ 0.3 AU. Based on these 
theoretically predicted orbital parameters, we discuss the detectability of binary planets by 
transit observations in the following sections. 


3. TRANSIT LIGHT CURVES OF BINARY PLANETS 

We discuss how the transit light curves of binary planets are modihed compared with 
those of a single planet, using the parameters given in the last section, with the aim of 
providing insight into the types of light curves that can be produced by a binary planet 
pair. For simplicity, we take the stellarcentric eccentricity of the binary barycenter as 
cg = 0 and the stellarcentric semi-major axis to be og = 0.4 AU, which is near the stability 
limit, OG.tide- Although transit detectability increases with decreasing og, the lifetime of 
binary planets is shorter. So, binary planets just outside aG,tide may be the most promising 
targets. 

We show the transit light curves of binary planets with equal-mass (Mp = Mj) and 
equal-radius (Rp). For the radius, we consider two cases: Rp = 2Rj and Rp = IRj. We 
set the orbital separation after the tidal circularization of binaries as au,i = 2.5Rtot, which 
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is a typical value obtained in Paper I. From Eq. o. the orbital separation after the orbit 
enters the spin-orbit synchronous state is estimated as au,2 — 3.9i?tot for Rp = 2i?j and 
Obi ,2 — 2.8i?tot for -Rp = 1-Rj. The corresponding binary rotation periods are ~ 5.4 and 3.3 
days for = 2i?j and R-p = li?j, respectively. Since Rq ~ 10i?j, the binary separations 
are ~ 0.8 and ~ 0.3 times the stellar diameter for Rp = 2i?j and Rp = i?j, respectively. 
As we pointed out in the last section, these binaries are most likely to be observed in the 
spin-orbit synchronous state. So, we assume a spin-orbit synchronous state for the binary 
planets. The stellarcentric orbital period is ~ 92 days, assuming a solar-mass central star 
and qq = 0.4 AU. For simplicity, we assume that the stellarcentric and binary orbits are 
coplanar. 


For light curves, limb-darkening is also taken into account, while stellar spots and 


pulsation are neglected. We calculate sample light curves u sing the method of 


using the quadratic limb darkening model of 


Pierce 


Pall (I2OI2I), 


fcOOOl) . For simplicity, we assume the 


radius, mass and luminosity of the host star are equal to the solar values of Rq, Lq, and 
Mq. 


Figure |2] shows an example of a transit light curve along with the positions of the 
binary planets, corresponding to the case where the binary’s barycenter passes the stellar 
center. In this case, the silhouette of the planets overlap (mutual transit) near the transit 
center. The top panel shows the transit light curve. The solid purple line represents the 
predicted light curve of a binary consisting of planets with Mp = IMj and Rp = 2i?j. For 
comparison, the light curve for a single planet with Mp = 2Mj, Rp = 2\/2Rj at semimajor 
axis a = 0.4 AU is also shown using a dotted red line. For this case we set Rp = 2\/2i?j, 
such that the transit depth is equal to the total transit depth of the two planets. Similarly, 
the light curve for a binary with Mp = IMj and radius Rp = Rj, and a single planet with 
2Mj and \/2i?j is also shown using a dashed blue line and a dash-dotted light-blue line 
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respectively. The middle panel gives the projected positions of the planets. The x-axis is 
the distance along the transit path with the origin at the center of the star. The bottom 
panel shows the projected positions of the transiting binary planets (the hlled black and 
gray circles) in the case of i?p = 2i?j for a number of snapshots during the transit. As 
shown in this illustration, the motion of the two planets around their common barycenter 
during the transit is taken into account when calculating the light curve. 

The leading planet of the binary (the hlled black circle) reaches the edge of the 
photosphere at f ~ —6 h from the transit center. Because of limb darkening and the area of 
the star being blocked by the planet increasing, the light curve gradually goes down with 
time. Because there is an off-set due to the binary separation, the ingress (the egress) is 
earlier (later) than for the single planet ht with a = 0.4 AU. For the case where only one 
transit is observed, and a single planet ht is done for the transit duration, the stellarcentric 
semimajor axis may be overestimated. The overestimation depends on the binary orbital 
phase during ingress and egress. Although overestimation is only slight in the case of Fig. [2l 
the maximum overestimation of the orbital period is a factor of 2, and thus, the semimajor 
axis would be overestimated by a factor of 2^/^ ~ 1.6. For the case where multiple transits 
are observed, the orbital period is known, and thus the semi-major axis can be accurately 
determined. 

At the transit center (f ~ 0), the light curve of the binary planet pair shows a bump 
due to a mutual eclipse between the planets. Depending on the positions of the two planets 
in their mutual orbit and orbital separation, a dip often appears at the transit center, 
instead of the bump (Figure |3]). Such a bump or a dip never occurs for a single planet 
transit, except for a bump resulting from a planet passing in front of a stellar spot. 

The most pronounced property of the light curves of binary planets is variation in the 
light curves from transit-to-transit. Figure |3] shows the light curves for a sequence of six 
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consecutive transits of the same binary system as that shown in Fig. [2j As can be seen, 
each light curve is different. The upper-left panel is the same as Figure |21 which we call 
“case A”. 

The three panels (the upper-right and upper/lower-middle panels) show a dip. Because 
the orbital separation is comparable to the stellar diameter for = 2i?j, the duration 
in which only one planet of the binary is transiting and that in which both planets are 
transiting are often comparable. In that case, the transit light curves show a deep dip 
(transit by both planets) near the transit center, sandwiched by relatively shallow transit 
(transit by one planet). We call this “case B”. In this case, the two-planet transit occurs 
without a mutual transit. However, since one or the other of the planets is transiting a 
highly limb darkened region, the transit depth is slightly shallower than for the single planet 
case. 

Also in the lower-right/left panels, the case of the transit of one planet and that of two 
planets can coexist. However, since the projected binary separations are smaller than those 
in case B, the curves show a “step” rather than a “dip.” Since the two planet transit is not 
affected by the limb darkening, the maximum transit depth is similar to the single planet 
£t. We call this case “case C”. 

Note that if the orbital separation is slightly larger, such that one planet enters the 
transit just at the time when the other leaves it, the transit curve has a bump at the transit 
center with a shallower transit depth than case A (Figure 0]). We call this case “case D”. 
We summarize cases A to D in table [TJ 

Note that if we consider the binary systems with non-equal i?p, the modulation of light 
curves becomes less pronounced. Figure Elis the same as Figure El except that the ratio of 
the physical radii of the binary planets is 2:1, keeping the total cross-section the same. In 
this non-equal case, statistical techniques for detecting exomoons may become necessary. 


13 


Case A is characterized by a mutual transit. From 
probability of binary planets with a mutual transit is 


Sato fc Asadal ( 2009 ). the detection 


P = P1P2, 


( 2 ) 


where pi = tobs/^K (^k is a stellarcentric Keplerian period of the binary barycenter) is the 
probability for the binary center to pass across the surface of the central star during the 
observational duration fobs, P2 = ^n/Tbi is the fraction of transits for which a mutual transit 
occurs during the transit duration = 2Rq/vk {P 2 = 1 for Tbi < fE), where Tb, is the 
Keplerian period of the binary system and uk is the stellarcentric Keplerian velocity of the 
binary barycen ter. For Uhi = 2.5(R i + Rj) = 10i?j, p ~ 0.1. So, the probability of case A is 


low. Note that 


Hirano et ah 


(1201211 detected a rare mutual transit, although it is for a pair 


of non-bounded planets. 


If the total mass of the binary system is determined through RV measurement, their 
mean bulk density can be calculated. However, as stated in the above, a single planet £t 
based on the transit duration leads to htting parameters of Mp = 2Mj and Rp = \/2i?j. 
As a result the calculated bulk density can be up to (\/2)^/2 = \/2 times smaller than 
the real one, with the exact factor depending on the derived impact parameter for both 
hts (see section SH] for an example). So, binary planets may be mis-classihed as inflated 
exoplanets, if a single planet £t is applied. The tidal stability limit ac.tide evaluated by the 
overestimated Rp is artihcially large, so if inflated exoplanets are located inside of ac.tide 
(but outside of oc.rui), h is worth considering the possibility that the object is really a 
binary planet pair. 


The changes in the light curves in a sequence of consecutive transits may be the most 
pronounced signal of planet binarity. However, if the stellarcentric orbital period is close to 
an integer number of times of the binary period and the total number of observed transit is 
relatively small, the transit-to-transit changes may not be signihcant. So, the detection of 
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perturbations to the light curve shape for each individual transit is also important. 

To see if such changes could be detectable, we have conducted a preliminary search of 
the open access CoRoT and Kepler transiting planet data. This investigation is discussed 
in the next section. 


4. DETECTABILITY OF BINARY PLANETS IN ARCHIVED COROT 

AND KEPLER DATA 


To investigate if binary planets analogous to those investigated in this paper could 
be detected using current technology, we focussed on data from the CoRoT and Kepler 
satellites. To provide context, both these missions, and the types of planets they were 
designed to detect will be discussed in turn. 


CoRoT is an ESA-led mission with the aim of using a 27cm diameter space telescope 
to detect transiting plan ets larger than Earth, in short period orbits as well as monitoring 


and characterising stars (lAuvergue et ah 


20091) . As the CoRoT satellite was in orbit of the 


Earth it suffered from thermal effects related to its orbital phase which led to predictable 
data errors. In addition, this satellite monitored a range of fields, located in the galactic 
centre and anti-center, with one observing run per held. Before the mission ended, nine long 
runs (150-90 days) and three short runs (20 days) towards helds in the galactic centre were 
completed along with six long runs, hve short runs, and one intermediate length run (50 
days) towards the galactic anti-center. As a result of this mission, ~ 100,000 light curves 
have been released and 27 planets have been detected. Given that less than a year’s worth 
of contiguous data was available for each candidate, it is unsurprising that most planets 
detected by CoRoT had small semi-major axes and short periods. 


On the other hand, Kepler is a Nasa-led mission with the aim of using a 0.95m 
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diameter telescope to discover an Earth-twin aronnd a Snn-like star flBorucki et al. 


2nn8h 


To ensnre less noise and a longer time baseline, Kepler was placed in an Earth-trailing 
orbit and monitored ~ 145,000 stars in one particular field. Apart from small gaps in 
the data due to satellite rotations, data downlinks and problems such as safe modes and 
coronal mass ejections, the data is continuous. As a result of the high data quality and the 
long time baseline, multiple transits for each candidate were routinely collected, allowing 
planetary orbital period to be derived and allowing sensitivity to sub-earth sized planets. 
These factors resulted in a much larger number of Kepler planet candidates (4234) as well 
as confirmed planets (978)1^ 


In this context, detection of binary planets using CoRoT and Kepler data will be 
discussed in turn. In particular, we show that simple light curve fitting programs can 
successfully identify binary planet candidates using a possible binary planet candidate from 
the CoRoT data set and confirmed binary star pair, transiting a brighter star, from the 
Kepler data set. Then, as the Kepler data set is a much richer place to search for binary 
planets, we conduct a simulation to demonstrate that Kepler quality data is sufficient to 
detect these planets. 


4.1. DETECTABILITY OF BINARY PLANETS IN COROT DATA 


To investigate binary planet detection, the CoRoT candidates were checked by eye. 


One inte resting object t 


detail in 


Erikson et al. 


rat w e discovered is CoRoT SRcOl E2 1066, which is described in 


(1201211 . The transit has a rel ative depth of 4 % an d duration of 66 


hours (see figure [6]) . From the long transit duration. 


Erikson et al. 


( 2012 ) suggested that 


this event might be a transit of an evolved or dwarf star by a distant gas giant planet. 


^Retrieved from http://kepler.nasa.gov/, 4*^ August, 2014. 
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where, by chance, the planet occulted a stellar active region (spot) at the centre of the 
transit. But this event could also be due to a transit of a binary gas giant planet pair with 
smaller stellarcentric semi-major axis. 


To check this claim we performed a single planet (dotted red) and binary planet (solid 
purple) fit to this data (see figure Ej). The single planet and binary planet models were 


calculated using the light curve simulation code of 


Pall fl2012!) . where, to account for long 


term trends in the light curve, the out-of-transit light curve was modelled using a cubic. 
For simplicity, we assume that the binary orbit normal is circular and is perpendicular to 
the line-of-sight and the chord made by the planet-planet barycenter across the star. The 
best fit models were then determined by locating the minimum residual using a simplex 
minimisation routine, with the fit improvement, calculated using the change in Bayesian 
Information Criterion (BIG), dehned as 


ABIC = nln 


/RSS 


- 


smg 


n 


In 


fRSSun 




n 


k ln(n) 


( 3 ) 


where RSSging and RSSbin are the residual sum of squares for the single planet and binary 
planet cases respectively, k is the number of extra degrees of freedom, five for this work, 
and n the number of data points used for the fit, in this case 1333. The derived parameters 
for the single and binary planet cases are shown in table [21 noting that the BIG value 
improves by ~ 80 when the binary, as opposed to the single planet model, is used, where 
a value of 10 is robust. As can be seen from table [21 the binary planet fit yields orbital 
and physical parameters for the putative binary tantalisingly similar to those derived in 
paper I. However, as only one transit is available (as this was a short run) it is not possible 
to rule out the starspot case. If more data were available for this target, its true nature 
could be determined by e.g. testing for long term trends due to light curve modulation due 
to starspots, resulting from stellar rotation. In addition, if another transit were detected, 
it would help differentiate between these cases by first, placing a strong constraint on the 
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orbital period (which is different in the binary and non-binary planet cases) and second 
indicating if there were transit-to-transit variations in transit width, depth and shape (see 
figures [3] to E]) as is likely for the binary planet case. These issues are discussed in greater 
detail in section 14.31 


4.2. DETECTABILITY OF BINARY PLANETS IN KEPLER DATA 


While no visually obvious binary planets are present in the K epler candidate and 


binary star data sets, a transiting binary star pair was discovered (jCarter et ah 


201111 . 


In this system, the central star has radius 2.0254 ± 0.0098 Rq while the transiting pair 
have radius 0.2543 ± 0.0014 Rq and 0.2318 ± 0.0013 Rq respectively, and the binary has 
a semi-major axis of 4.729 Rq. As the host star is between 3000 and 5000 times brighter 
than the members of the transiting binary pair, and the radius ratios are similar to those of 
the planetary binaries described in this paper and in paper I, these light curves should be 
approximately analogous to those of the planetary case, and we can test our fitting code on 
a real transiting binary system. 


Example transits from this system are shown in figure [71 Note how they look markedly 
different from those from a single object transit and similar to some of the sample transits 
shown in figures [3] and jH We fit this data with our simplified fitting model, given by the 
solid (purple) line. As this transiting binary system has an inclined orbit, and the impact 
parameter for the transit is large, we do not expect the light curve and the model to exactly 
match. In particular, this binary orbit is inclined such that the leading star transits closer 
to the stellar limb than the trailing star and is the reason why the transit duration and 
depth observed for the leading star is shorter and shallower than predicted by the fit and 
the transit duration and depth observed for the trailing star is longer and deeper than 
predicted by the fit. While the fit isn’t perfect, the binary system is unequivocally detected 
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with a reduction in the BIG of over 2000, which indicates that our simple model still can 
successfully detect binary systems for inclined binaries and high transit impact parameter 
systems. 

To investigate if binary planets analogous to those investigated in this paper could be 
detected in real Kepler observations, a preliminary investigation was conducted. Focussing 
on the long cadence data, we selected a Kepler candidate, KOI 3681.01 (KIC 2581316), 
showing an anomalously large radius (22 Earth radii), transit duration (21.3 hr) and with a 
host magnitude close to the 12^^ magnitude Kepler target specifications (11.69), with the 
aim of demonstrating that it is possible to robustly detect the difference between a binary 
planet transit light curve and a transit light curve due to an inflated single planet. To do 
this, we simulated realistic noisy binary planet transit light curves using two different noise 
sources, KIC 2581316/KOI 3681, a 11.69 magnitude star and KIC 9517393/KOI 2076, a 
15.3 magnitude host star, and then tried to detect the planet binary by fitting a single 
planet and binary planet light curve model. The method used to produce these light curves 
is described in the next section. 


4.2.1. Simulating Noisy Binary Planet Lightcurves 


Simulated realistic binary planet light curves were constructed from real Kepler light 
curves in five stages. First, long cadence data was downloaded from the NASA Exoplanet 
Archivqj corresponding to quarters 0 to 16. Second, all planetary transits were completely 
removed. Third a sequence of binary planet transits was simulated. Fourth, out-of-transit 
data was used to give the simulated transits realistic noise. Finally, the simulations and 
the original data were stitched together, where a cubic was added to the simulation to 


“^http://exoplanetarchive.ipac.caltech.edu/index.html 
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ensure that the gradient and value of the endpoints matched. The hnal four stages will be 
discussed in detail. 


For both KIC 2581316 and KIC 9517393, all planetary transits were first identified. To 
ensure that no signal corresponding to real moons of planetary candidates in these systems 
remained, estimates of the planetary masses were used to calculate reasonable values for 
the Hill sphere, and all data corresponding to the transit of the planet or the Hill sphere 
was removed. 


Then, using the code of 


Pall (120121 ). we simulated sequences of transits of binary 


planet pairs. Following Paper I, we simulated two classes of binary planet pairs, one where 
both components had equal radii and one where the radii of one of the components was 
twice that of the second. To ensure that the transit depth for the binary case was similar 
to that observed for the real candidate (see £gure[8]), for the equal radius case we set 
Ri = R 2 = 0.06297?* while for the other case we set Ri = 0.07957?* and 7?2 = 0.03987?*. 
Also, the semi-major axis of the planet binary was taken to be 7?* corresponding to a 
typical binary planet separation (see section [3]), and the transit velocity was altered to 
approximately match the transit duration. In addition, following the analysis in section [3l 
the planets are assumed to be equal mass and, for simplicity, these systems are assumed to 
have zero eccentricity, be coplanar with the planet orbit and have impact parameter equal 


to that of KOI 3681.01. 


To ensure these light curves displayed realistic photometric noise, we randomly selected 
sections of out of transit light curve from Ql-15 data from KIC 2581316 for the case of low 
noise and KIC 9517393 for the case of high noise. Our simulated light curves were then 
multiplied by these sections of of data to produce noisy binary planet light curves. 

Finally, the simulated light curves were stitched into the original data. To ensure that 
the light curves were continuous, a cubic was added to each transit light curve such that the 
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gradient and value at both edges of the simulated transit light curve matched the gradient 
and value at the edges of the original data. 


Using this method we simulated 50 six transito sequences of light curves for these 
systems (see figure [8] for examples), where the initial true anomaly was randomised for each 
sequence. 


4-2.2. Determining if the Binary Planets were Detectable 

To determine if planet binarity was detectable in these simulated light curves, using 
our light curve fitting code, we fitted each sequence of transits with a single planet and 
binary planet model and recorded the difference in the BIG. To determine the effect of data 
length we repeated the process including only one and only three transits. Some example 
fits are shown in figure [H] and the results are shown in figures [9] and [TOl All binary planets 
were detected, nearly all, robustly. 

To investigate the behaviour of the detection threshold as a function of host star 
magnitude, three additional stars were chosen from the Kepler catalog, KIC 12121701 
(magnitude 15.611), KIC 8827930 (magnitude 15.999) and KIC 2438406 (magnitude 
16.546). While these results apply to these specihc systems, the trend should be indicative 
of the true trend. The analysis described previously was repeated for the equal radius ratio 
and one or three transit cases, with the results plotted in hgure [TT] along with the results 
for KIC 9517393. As can be seen, in our sample, binary planets are robustly detected for 
stars with magnitude less than 16 and 15.5 for the three and one transit cases. 

In addition, binary planets may exist with radius ratios outside of those investigated in 


^KOI 3681.01 shows seven transits while KOI 2076.02 shows six. 
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Paper I. To investigate the detectability of such systems, simulations with a range of radius 
ratios were constructed using photometric noise from KIC9517393, the 15.3 magnitude 
target, and analysed using our code. As can be seen from £gure[T21 for the case where one 
transit is observed, binary planets are robustly detected for radius ratios smaller than 2:1 
while for the case of three observed transits, even planet pairs with radius ratio 5:1 are 
robustly detected. 

From these simulation results we show that a range of binary planets, analogous to 
those simulated in this paper could be practically detected in Kepler data even for the case 
where the star is dim (15^^ magnitude) and the number of transits is small. In particular 
we show that binary planets are more detectable around host stars with low relative 
photometric noise compared to high photometric noise targets. In addition, we show that 
detectability improves when the number of observed transits increases. 

Finally while this analysis does include many important physical factors e.g. realistic 
photometric noise, transit to transit variation, it does not include the effect of spot crossing 
events in the light curve. As discussed in sectional possible spot crossing events may closely 
mimic the effect of gas giant binary planet transit, so this topic warrants discussion. 


4.3. THE EFFECT OF STARSPOTS ON BINARY PLANET DETECTION 

One factor which may hamper the detection of binary planets is the presence of spot 
crossing events in a transit light curve. We propose that it should be possible to differentiate 
between a single planet which orbits a spotty star and a binary planet given sufficient 
numbers of transits by either conhrming the presence of spots or by providing supporting 
evidence for a binary planet. In particular this can be done in a number of ways including: 

1. Investigate if the timing of spot crossing events corresponds to a physically realistic 
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star: For the case where the planetary orbit is inclined with respect to the stellar 
spin axis, the transit chord may cross one or more active latitndes. As a resnlt, spot 
crossing events are more likely to appear on the same part of different transits light 
cnrves, with the position corr esponding to the location of the active latitndes (e.g. 


Sanchis-Ojeda fc Winn 


20111 ). In addition, it has been snggested that interactions 


between the planet and the star can also lead to correlation between the position 


of the planet and the position of act ive regions. T 


and detecta ble in some systems (e.g. 


others (e.g. 


Miller et al 


20121 : 


Pagano et al 


Scandariato et al. 


lis phenomenon may be 


2009 


Herrero et al. 


present 


2013 ) bnt not 


20131) . however, as the planets of 


interest for this work are distant (ontside 0.3AU), snch an interaction is likely to be 
very weak or absent. Conseqnently, the presence of light cnrve pertnrbations that 
cannot be predicted, bnt occnr with higher probability in certain sections of the light 
cnrve, conld mean strong evidence for starspots being the canse. However for the case 
of binary planet pairs, the size and location of bnmps in the light cnrve shonld relate 
to the sizes of the planets and will always correspond to a physically realistic planet 
model. 


2. Use trends in the out-of-transit light curve to co 
canse changes to the ont-of-transit light cnrve. 


nfirm presence of spots: Sta rspots also 


Sanchis-Oieda et al. 


( 20121) snggested 


corellating long term light cnrve modnlations and spot crossing events to determine 
relative inclination in mnlti planet systems, bnt the same processes conld be nsed 
to provide evidence for or against a spot being the canse of a particnlar light cnrve 
featnre. 


3. Image the stellar surface to confirm presence of spots: For a possible evolved host 
star, they are known to host large starspots. Snch spot s have been imaged (e.g., 

3 


Vogt fc Penroc 


19§3; 

Vogt et al. 

19§7; 

Strassmeier 


2002! ). 


For some host stars this 
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may be an option to determine the presence of spots. 

4. Compare measured planetary mass to predicted mass: As mentioned previously, for 
the single planet case, the mass derived from radial velocity is likely to match the 
measured radii, while for binary planets that have been incorrectly classihed, it will 
be systematically low. 

Similar to the simulations presented, these arguments indicate that the detectability of 
binary planet pairs increases as the number of transits increases. In addition, observations 
that are likely to be taken to try to prove planetary nature e.g. RV observations can also 
be used to determine likelihood of a given candidate being a binary planet. 


5. CONCLUSIONS 

In this paper, we have discussed the possibility of the observational detection of 
extrasolar binary planets (gravitationally bound pairs of gas giant planets) by transit 
observations, based on the results of N-body simulations on the tidal capture between gas 
giants and calculations of long-term tidal evolution after the capture performed by Paper I. 
Paper I showed that the formation probability of a planetary binary is as much as ~ 10% 
almost independent of stellarcentric semimajor axis, of the binary (og) and predicted that 
the typical binary separation is 3-5 times the sum of physical radii of the planets and the 
binary planets are tidally stable for 10^° years if og > 0.3 AU. 

Using these constraints, we have modelled transit light curves of physically plausible 
binary planets. These light curves have a deep dip, a big bump or a step or two separate dips, 
and are noticeably and statistically different from those of a single planet. Furthermore, the 
transit shape changes from transit to transit, compared to the single planet case. Because 
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of these features, the transits of binary planets might be classihed as false positivesjj If RV 
measurement is also available, the bulk density can be estimated. The single planet £t for 
two equi-sized binary planets can give a lower bulk density than the real value. Thus true 
binary planets could also have been classihed as “inhated” planets if oq > 0.3 AU. 


We show that the CoRoT target SRcOl E2 1066 is well £t by a binary planet model 
and put forward the alternate scenario that it could als o be due to th e tran sit of a binary 


planet in addition to the starspot scenario proposed by 


Erikson et ah 


1 20121) . In addition. 


we show that binary planets may be present in, and would be detectable in the Kepler data 
set and are most detectable where the host star shows little noise and a number of transits 
are available. In addition, for host stars with magnitude less that 15, we show that a broad 
range of binary planets are robustly detectable, even for the case of one or a few observed 
transits. Prompted by this preliminary analysis we propose that we should do an accurate 
reanalysis of the irregular changing light curves of Kepler and CoRoT planets and planet 
candidates orbiting beyond 0.4AU from their central star, where over 100 candidates and 
false positives exist. 


We are grateful for the offer of the gravitational lensing calculation program by 
Takahiro Sumi. We also thank Tristan Guillot and Rosemary Mardling for discussions on 
observations of binary planets. In addition, we thank Jessie Christiansen for advice on the 
Kepler pipeline. We would also like to thank an anonymous referee for helpful comments 
which improved the quality of this paper. This research was supported by a grant for JSPS 
(23103005) Grant-in-aid for Scientihc Research on Innovative Areas. KML was supported 


Light curves that change irregularly have been already discovered (e.g.. 


Barnes et ah 


201311 . although this object is inside of the tidal stability limit and would not be a binary 
planet. 
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Fig. 1.— The distributions of binary separations. The left panel shows the separations (obi) 
just after tidal trapping followed by binary orbital circularization by planet-planet dynamical 
tides that are obtained from the N-body simulations conducted in Paper I. The separations 
(abi,i) after entering the spin-orbit synchronous state as a result of long-term quasi-static 
tidal evolution are plotted in the right panel. 
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Fig. 2.— The top panel shows transit light cnrves for simnlated binary planets transiting a 
solar-like limb darkened star for the case where the binary’s barycenter passes the center of 
the photosphere as the binary planet nndergoes a mntnal event along the line-of-sight. The 
cases of a IMj, 2i?j binary (solid pnrple line), 2Mj, 2\/2i?j single planet (dotted red line), 
IMj, IRj binary (dashed bine line) and 2Mj, \/2i?j single planet (dash-dotted light-bine 
line) are shown. The middle panel gives the positions of the members of the 2i?j (solid 
pnrple line) and Rj (dashed bine line) binary pairs along the transit chord as a fnnction 
of time. The bottom hgnre shows the positions of the transiting binary planet (the hlled 
black and grey circles) in front of the photosphere (the hlled yellow circle). The x-axis is in 
the direction of motion of the binary’s barycenter and normalized by the Snn’s radins, Rq. 
The origin of x-axis is taken at the center of photosphere. The barycenter of binary always 
remains on the x-axis and the two planets move in coplanar orbits. 
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Fig. 3.— A sequence of six consecutive transit light curves and the projected positions of 
the transiting binary planets, for the case of Fig. |2l The upper-left panel is the same as 
Fig. m We set t = 0 when the binary’s barycenter passes the stellar surface center in all 
hgures. The light curve in the upper-left panel is case A, the light curves in the upper-right 
and upper/lower-middle panels are case B, and those in the lower-right/left panels are case 
C. 
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Fig. 4. — An example of a transit light cnrve of a binary planet in case D. We set = 2i?j 
and the orbital separation that is 1.25 times larger than that in Fignre[3](~ 3.9i?tot)- 
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Fig. 5.— Same as Figure [HI except that the ratio of the physical radii of the binary planets 
is 2:1, keeping the total cross-section the same. 
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Fig. 6.— Light curve for CoRoT target SRcOl E2 1066, where only the data points (black 
crosses) with status flag set to 0 are plotted. The single planet (dotted red) and binary 
planet ht (solid purple) are also shown. The best ht values corresponding to these fits are 
given in table ID 


case 

transit depth 

transit duration 

transit shape 

A 

half 

similar 

bump due to mutual transit 

B 

slightly shallower 

longer 

deep dip 

C 

similar 

longer 

asymmetric step 

D 

half 

twice 

bump due to wide separation 


Table 1: The light curve properties of a binary planet pair compared with those of an anal¬ 
ogous single planet with the same stellarcentric semi-major axis and total cross-sectional 


area. 
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Single planet fit 


Impact parameter 

0.573043 

Transit velocity (i?*/day) 

0.738503 

Semi-major axis (AU) 

25 

Planet radius (i?*) 

0.187059 

Binary planet fit 

Impact parameter 

0.847337 

Transit velocity (i?*/day) 

1.044480 

Semi-major axis (AU) 

12.5 

Planet 1 radius (i?*) 

0.22087 

Planet 2 radius (i?*) 

0.156546 

Binary semi-major axis {R^,) 

0.984286 

Mid-transit true anomaly 

1.216958 


Table 2: The best fit values corresponding to the single planet and binary planet hts shown in 
£gure[6l The semi-major axis listed is calculated assuming a Sun-like h ost star and a c ircula r 
orbit of the planet or planet binary around the star. If, as suggested by 


Erikson et ah 


( 120121 ), 


the host star is evolved, the semi-major axes calculated for the fits will be much lower. 
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Time (days) 


Fig. 7.— The first six transits of the binary star pair in the KOI-126 system (black crosses) 
along with the best £t single star (dotted red line) and binary star (purple line) model 
calculated using our simplified fitting code. 
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Fig. 8.— The real transits of KOI 3681.01 (left) from quarters 2, 5, 7, 9, 12 and 14, along 
with simulated transits for the case of a binary pair with radius ratio 2:1, an equal radius 
binary pair, and an equal radius binary pair where out of transit data from KIC 9517393 
was used to contaminate the light curve (right). The single planet (dotted red) and binary 
planet (solid purple) hts are shown for the simulated data. 
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A BIC 



A BIC 


Fig. 9.— Distribution of the decrease in BIC for the set of simulations of an equal radius 


binary planet along with the robust detection limit 
(red dashed line). 





































Frequency Frequency 


39 



A BIC 



Fig. 10.— Distribution of the decrease BIC for the set of simulations of a binary planet with 
radius ratio 2:1, again plotted with the robust detection limit (red dashed line). 
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Fig. 11.— Box plot showing the minimum, hrst quartile, median, third quartile and maxi¬ 
mum decrease in BIC for the set of simulations of an equal radius binary planet for a range 
of different stellar magnitudes along with the robust detection limit (red dashed line). 
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Fig. 12.— Box plot showing the minimum, hrst quartile, median, third quartile and max¬ 
imum decrease in BIC for the set of simulations of binary planets with a range of radius 
ratios with noise taken from KIC 9517393 with the robust detection limit (red dashed line). 














































